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We obtain the constraints on the diffuse neutrino background from primordial black hole evapo- 
rations using the recent data of KamLAND experiment setting an upper limit on the antineutrino 
flux from the Sun. 



I. INITIAL PBH MASS SPECTRUM 

' According to Press-Schechter model, the fraction of mass in the universe contained in gravitationally bounded 
' object more massive than M is [2] 
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. Here, j^dM is the number density of objects of mass in (M, M + dM), and the integral on the left-hand side gives 
the total mass density in overdense regions more massive than M. P{5R,M,ti) is the probability that a randomly 
^-H ■ chosen point in the universe has density contrast 6r > 8c at time U, and M is the smoothing scale of 5r {R is the 
^ ' radius of the spherical overdense region) . The value 5c is some critical value of the density contrast which corresponds 
, to the beginning of the non-linear regime. In a case of a gaussian field of density fluctuations this probability is 

^ : P- f-^) , (1.2) 



] where = a\{M,t) —< Sj^ > is the variance of the smoothed density. The mass distribution of gravitationally 

■ bounded objects is [2] 

Q ■ dM M dM' ^ ' ' 

' The variance ct|j is given by the relation [3] 
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where W{kR) is the smoothing window function, Snik,!) is the horizon crossing amplitude, T(fc) is the transfer 
function. 

For the calculation of the mass spectrum of PBHs one must use the connection between the PBH mass Mbh and 
characteristics of the overdense region. In the model of the near-critical collapse proposed in ref. [11] one has the 
relation: 

Mbh = M^/'M2/3fc(<5g - ScV", (1.5) 

where Mi is the horizon mass at time ti when the density fluctuations develop, and 6^ is the smoothed density 
contrast at the moment when the fluctuation crosses horizon, 

6^ = {M/M,f/'6R. (1.6) 

In the approach of ref. [4] all PBHs have masses roughly equal to the horizon mass at a moment of the formation, 
independently of the perturbation size. The approximate connection between the PBH mass and the horizon mass in 
such a model is very simple: 

Mbh = 7'/'M/^'Af2/3. (1.7) 
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where 7 is the ratio of the pressure to energy density (7 = 1/3 in the radiation dominated era). 
Using these formulas one obtains the final expression for the PBH mass spectrum [6] 
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Here we use the general notation for the connection between Mbh, M and i5^, 

Msff = /(M,5i). 

Particular cases of this connection are given by Eqs.(1.5,1.7). 
Substituting the expression for the window function, 



W{kR) = 3 



/ sin kR cos kR 
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and for the transfer function [5], 



and the parametrization 



SH{k) = Sniko) - 
\ko 

for the horizon crossing amplitude into Eq.(1.4) for the variance cr^ one obtains, approximately. 
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The exact value of the coefficient depends on the spectral index (for n = 1.3 = 6.82). In Eq.(1.13) kfi is the 
comoving wave number, characterizing the perturbed region, kfi = ^, and M/, is the fluctuation mass at the moment 
when the perturbed region enters the horizon at the radiation dominated era. The value of is simply connected 
with kfi: 



kfi = tteqllgq 
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At the same time, an initial value of the fluctuation mass, M, is proportional to R^, so Mh ^ M^/^. The exact 
connection is 



The denominator {aHY in Eq.(1.13) must be taken &tt = ti, and one can easily see that 
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At last, in the spectrum formula (1.8) the following notation is used 
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The limits of integration in Eq.(1.8) are determined using the conditions 

S^^^n=S, , Mmin=Mr^ =Mi . (1.18) 

We have, for calculations oinBH^MBH) two free parameters: the spectral index n, giving the perturbation amplitude 
through the normalization on COBE data on large scales, and ti, the moment of time just after reheating, from which 
the process of PBH formation started. The value of ti is connected, in our approach, with a value of the reheating 
temperature, 

t.=0.301,.-^/^|^ (1.19) 
(5* 100 is a number of degrees of freedom in the early Universe). The initial horizon mass M, is expressed through ti: 

M^ = ^^Tt^,P^■ (1.20) 

In the case of Page-Hawking [4] collapse the minimum mass in the PBH mass spectrum is given by the simple 
formula 

M™™ = (1.21) 

which follows from Eq.(1.7). Correspondingly, the minimum value of PBH mass is determined entirely by the reheating 
temperature. In the case of the near-critical collapse there is no minimum value of the PBH mass in the spectrum 
because 5^ can be almost equal to (see Eq.(1.5)). In this case the PBH mass spectrum has a broad maximum at 
Mbh ~ Mi. 



II. NEUTRINO DIFFUSE BACKGROUND FROM PBHS 



Evolution of a PBH mass spectrum due to the evaporation leads to the approximate expression for this spectrum 
at any moment of time: 

2 

"^^^""'^^ = (3a^ + m3)2/3 "^^ ((3ai + m3)V3) , (2.1) 

where a accounts for the degrees of freedom of evaporated particles and, strictly speaking, is a function of a running 
value of the PBH mass m. In all our numerical calculations we use the approximation 

a = const = a(Mg'5^), (2.2) 

where M^^^ is the value of Mbh in the initial mass spectrum corresponding to a maximum of this spectrum or, 
in the case of Page-Hawking collapse, to a minimum value of PBH mass in the spectrum. Special study shows that 
errors connected with such an approximation arc rather small. The detailed calculations of the function a{m) were 
carried out in the works [14]. Here we use the simplified approach in which a{m) is represented by a series of step 
functions [8] 

The expression for a spectrum of the background radiation is [8] 
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(2.3) 



In this formula a^, a and aq arc cosmic scale factors at i^, t and at present time, respectively, and f{E,m) is a 
total instantaneous spectrum of the background radiation (neutrinos or photons) from the black hole evaporation. 
It includes the pure Hawking term and contributions from fragmentations of evaporated quarks and from decays of 
pions and muons (see [8] and earlier papers [14-16] for details). 

In the last line of Eq.(2.3) we changed the variable t on z using the flat model with JIa = for which 

|=-^^(^.™(.+i)3+..(.+i)r^^ 

(2.4) 

nr = (2.4 • IQt^h?)-^ , h = 0.67. 
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III. RESULTS 



In some theoretical schemes (e.g., in the model of spin-flavor precession in a magnetic field) the Sun can emit a large 
flux of antincutrinos. The LSD experiment [18] sets the upper limit on this flux, < 1.7 %. The corresponding 

constraints on diffuse neutrino background from primordial black holes were obtained in the work of authors [8]. 

Recently, the KamLAND experiment obtained a new limit on electron antineutrino flux from the Sun, according 
to which this flux is less than 2.8 x 10~^ of the SSM v^. flux. The neutrino detection was carried out in both 
experiments using the reaction 

j>e +p ^ n + e+. (3.1) 

No candidates were found in [1] for an expected background of 1.1 ± 0.4 events. We calculate the limit on the 
antineutrino flux using the formula 

= ^■A^.y.Ar,.4. ' ^^-2) 

where a is the cross-section of the reaction (3.1) averaged over neutrino background spectrum from PBHs, Ai? is the 
neutrino energy range measured in the experiment (8.3 — 14.8 MeV), T = 1.6- 10^ s is the livetime, and ISS^ = 4.6- 10^^ 
is the number of target protons in the fiducial volume. 

Using the results of calculations of electron neutrino background spectra from PBHs (the examples of such spectra 
are given on Fig.l for the case of Page-Hawking collapse) and normalization of the perturbation amplitude SH{k) on 
COBE data (for details see [8]), we obtain the constraints on the spectral index n shown on Fig. 2. 
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FIG. 1. Electron neutrino background spectra from PBHs, calculating for several values of the spectral index. Dashed curve 
shows the theoretical atmospheric neutrino spectrum at Kamiokande site [17] (averaged over all direction) 
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FIG. 2. Constraints on spectral index from KamLAND experiment data (solid lines) in comparison with constraints from 
data on cxtragalactic gamma ray background (dotted lines). Upper lines correspond to critical collapse scenario of the PBH 
formation, lower ones - to standard case. 
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